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In  this  paper  we  investigate  the  growth  rates  of  Gortler  vortices  in  a  compressible  flow 
in  the  inviscid  limit  of  large  Gortler  number.  Numerical  solutions  are  obtained  for  0(1) 
wavenumbers.  The  further  limits  of  (i)  large  Mach  number  and  (ii)  large  wavenumber  with 
0(1)  Mach  number  are  considered.  We  show  that  two  different  types  of  disturbance  modes  can 
appear  in  this  problem.  The  first  is  a  wall  layer  mode,  so  named  as  it  has  its  eigenfunctions 
trapped  in  a  thin  layer  near  the  wall.  The  other  mode  we  investigate  is  confined  to  a  thin 
layer  away  from  the  wall  and  termed  a  trapped  layer  mode  for  large  wavenumbers  and  an 
adjustment  layer  mode  for  large  Mach  numbers,  since  then  this  mode  has  its  eigenfunctions 
concentrated  in  the  temperature  adjustment  layer.  We  are  able  to  investigate  the  near 
crossing  of  the  modes  which  occurs  in  each  of  the  limits  mentioned. 
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tract  No.  NASl-18605  wliile  the  authors  were  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23(365. 
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1.  Introduction 


Our  aim  in  this  paper  is  to  investigate  the  growth  rates  of  Gortler  vortices  for  a  com¬ 
pressible  inviscid  flow  over  an  infinite  cylinder  in  the  limits  of  high  Mach  number  and  high 
wavenumber.  This  investigation  is  motivated  by  recent  interest  in  the  development  of  hyper¬ 
sonic  aircraft  which  might  well  be  capable  of  reaching  speeds  in  the  order  of  20  -  25  Mach. 
Real  gas  effects  will  certainly  come  into  play  at  these  speeds  but  for  simplicity’s  sake  we  have 
not  taken  them  into  account  in  this  paper.  We  have  also  assumed  that  Chapman’s  viscosity 
law  holds  for  this  fluid. 

The  most  obvious  difference  between  Gortler  vortices  in  incompressible  and  hypersonic 
flows  is  that  the  presence  of  a  temperature  adjustment  layer,  where  the  temperature  decays 
rapidly  to  its  free  stream  value,  at  the  edge  of  the  boundary  layer  enables  hypersonic  Gortler 
vortices  to  be  concentrated  well  away  from  the  wall.  In  the  incompressible  case,  we  know, 
from  the  work  of  Hall  (1982a, b,  1983),  and  Denier,  Hall,  and  Seddougui  (1991),  that  unstable 
Gortler  vortices  are  not  localized  within  the  boundary  layer  for  order  one  Gortler  numbers. 
For  higher  Gortler  numbers  the  most  dangerous  Gortler  vortices  have  wavelengths  small 
compared  to  the  boundary  layer  thickness  and  are  trapped  near  the  wall.  Not  surprisingly 
this  situation  does  not  change  significantly  for  order  one  Mach  numbers  and  this  nonparallel 
problem  has  been  discussed  by  Wadey  (1990)  and  Spali  and  Malik  (1989).  In  the  latter 
two  papers,  the  nonparallel  equations  were  solved  numerically  following  the  approach  of  Hall 
(1983)  and  the  main  result  obtained  was  that  the  growth  rate  of  a  Gortler  vortex  is  a  function 
of  its  upstream  history.  The  numerical  calculations  of  Wadey  (1990)  do  suggest  that  as  the 
Mach  number  increases  the  unstable  Gortler  vortices  locate  themselves  towards  the  edge  of 
the  boundary  layer.  This  view  is  supported  by  Hall  and  Fu  (1989)  whose  main  result  was  that 
the  logarithmically  small  temperature  adjustment  layer  at  the  edge  of  a  hypersonic  boundary 
layer  can  support  Gortler  vortices  and  the  most  dangerous  wavelengths  of  the  vortices  are 
comparable  with  the  thickness  of  this  layer.  Fu,  Hall,  and  Blackaby  (1990)  have  considered 
the  influence  of  real  gas  effects  and  Sutherland’s  viscosity  law  on  the  Gortler  instability  in 
hypersonic  flows. 

This  present  paper  is  restricted  to  linear  regimes  of  vortex  growth,  for  a  detailed  account 
of  nonlinear  regimes  the  reader  is  referred  to  the  review  article  by  Hall  (1990).  It  is  also 
worth  noting  that  hypersonic  boundary  layers  are  susceptible  to  instabilities  other  than 
those  caused  by  streamline  curvature,  such  eis  ToIImien-Schlichting  wave  instabilities  which 
have  been  discussed  by  Cowley  and  Hall  (1990)  and  Smith  and  Brown  (1990).  Clearly  any 
nonlinear  investigation  of  Gortler  vortices  at  hypersonic  speeds  must  allow  for  the  possible 
interaction  of  the  vortices  and  these  other  types  of  instabilities. 

The  present  paper  is  concerned  with  the  inviscid  limit  of  Gortler  vortices  in  a  compressible 
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boundary  layer.  We  find  that  two  distinct  modes  exist  and  we  are  able  to  describe  them 
in  the  separate  limits  of  large  free-stream  Mach  number  and  large  disturbance  wavenumber 
(although  for  the  case  of  large  free-stream  Mach  number  the  disturbance  wavenumber  will 
also  be  large).  The  first  mode  also  exists  in  an  incompressible  boundary  layer  and  is  termed 
a  wall  layer  mode  here  since  its  eigenfunctions  are  concentrated  in  a  thin  laj^er  near  to  the 
wall.  For  the  incomprecsible  case,  in  the  limit  of  large  wavenumber,  this  mode  is  described  by 
Denier,  Hall,  and  Seddougui  (1991).  The  second  mode  has  its  eigenfunctions  concentrated  in 
a  thin  layer  away  from  the  wall  and  therefore  is  referred  to  as  a  trapped  layer  mode  in  the  large 
wavenumber  limit.  In  the  hypersonic  limit  this  layer  is  precisely  the  temperature  adjustment 
layer,  a  logarithmic  layer  where  the  basic  temperature  changes  ranidlv  from  its  0(  Dj^)  value 
close  to  the  wall  to  its  0(1)  value  at  the  edge  of  the  boundary  layer.  Consequently,  in  the 
hypersonic  limit  this  mode  is  termed  an  adjustment  layer  mode. 

The  numerical  solution  of  the  equations  governing  the  stability  of  Gdrtler  vortices  in  the 
inviscid  limit,  discussed  in  Section  2,  show  that  the  growth  rates  of  the  infinity  of  solutions 
of  the  two  modes  described  above,  as  functions  of  the  disturbance  wavenumber,  appear 
to  intersect.  A  similar  near-crossing  of  modes  is  evident  in  the  numerical  results  of  Mack 
(1987)  in  the  form  of  kinks  in  the  neutral  curves  of  the  generalized  inflection  point  modes 
for  compressible  flow  over  a  flat  plate.  Asymptotic  solutions  for  the  infinity  of  solutions  of 
the  compressible  Rayleigh  equation,  termed  acoustic  modes,  in  the  hypersonic  limit  have 
been  given  by  Cowley  and  Hall  (1990),  while  the  so-called,  single,  vorticity  mode  for  large 
values  of  the  Mach  number  has  been  discussed  by  Smith  and  Brown  (1990).  Cowley  and 
Hall  (1990)  postulated  that  the  near-crossing  of  the  neutral  curves  could  be  described  by 
a  WKB  description  of  the  acoustic  modes  and  the  vorticity  mode.  This  is  precisely  the 
method  employed  by  Smith  and  Brown  (1990)  in  their  investigation  of  inviscid  modes  of 
instability  for  large  Mach  number  flows.  The  results  of  Smith  and  Brown  (1990)  show  that 
the  discontinuous  vorticity  mode  becomes  continuous  in  the  limit  of  large  Mach  number. 
Additionally,  they  ascertain  that  for  large  Mach  numbers  the  acoustic  and  vorticity  modes 
are  separated  by  an  exponentially  small  amount  as  was  proposed  by  Cowley  and  Hall  (1990). 

The  objective  of  the  present  paper  is  to  describe  the  wall  layer  modes  and  trapped  layer 
modes  present  for  inviscid  Gdrtler  vortices  and  to  investigate  their  near-crossing  which  exists 
as  outlined  above.  To  this  aim  we  follow  the  ideas  of  Cowley  and  Hall  (1990)  and  Smith  and 
Brown  (1990)  and  consider  a  WKB  analysis  of  the  modes  of  interest. 

The  layout  of  this  paper  is  as  follows.  In  Section  2  we  derive  the  equation  governing  the 
structure  of  inviscid  Gdrtler  vortices  in  a  compressible  boundary  layer  and  then  discuss  some 
numerical  results  of  this  equation.  In  Section  3  we  take  a  closer  look  at  the  wall  layer  modes, 
those  with  their  eigenfunctions  trapped  near  the  wall,  in  thu  hypersonic  limit.  We  follow 
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this  by  looking  at  the  adjustment  layer  modes,  those  whose  eigenfunctions  are  concentrated 
in  the  temperature  adjustment  layer,  in  the  hypersonic  limit  in  Section  4.  Then  in  Section 
5  we  look  at  the  wall  layer  modes  for  the  high  wavenumber  limit.  In  Section  6  we  consider 
the  trapped  layer  modes  for  the  high  wavenumber  limit,  these  are  the  equivalent  of  the 
adjustment  layer  modes  for  the  hypersonic  limit.  We  then  consider  the  near  crossing  of  the 
different  types  of  modes  for  both  limits  in  Section  7  before  finally  summarizing  our  results 
in  Section  8. 

2.  Formulation 

Our  aim  in  this  section  is  to  obtain  the  equation  which  determines  the  structure  of  a 
Gortler  vortex  in  a  compressible  boundary  layer.  The  boundary  layer  considered  is  that  of 
a  flow  over  the  cylinder  y*  —  0.  —oo  <  z*  <  oo  so  that  the  2*-axis  is  a  generator  of  the 
cylinder  and  y*  measures  the  distance  normal  to  the  surface.  The  a:’-coordinate  measures 
distance  along  the  curved  surface,  which  is  supposed  to  have  variable  curvature  j J 1) 
where  m  and  I  are  length  scales.  The  Reynolds  number  R,  Gortler  number  G  and  curvature 
parameter  8*  are  defined  by 


V 

(2.1a) 

G  =  2Rh\ 

(2.16) 

5*  =  — , 

m 

(2.1c) 

where  Uqo  is  a  typical  flow  velocity  in  the  streamwise  direction  and  v  is  the  kinematic  viscosity 
of  the  fluid.  The  Reynolds  number  is  assumed  to  be  large,  whilst  8*  is  sufficiently  small  so 
that  as  >  0  the  parameter  G  is  fixed  and  of  order  one.  We  take  the  basic  two-dimensional 
boundary  layer  to  be  of  the  form 


u=:Uoo[u{X,Y),R-h{X,Y),0][l 

+  o(/?-h], 

(2.2) 

where 

X*  y*R^ 

X-  Y-  ^  . 

(2.3) 

We  chose  to  look  at 

a  Blasius  boundary  layer  by  putting 

tZ  =  /'(//)  where  /  satisfies 

2/'"  +  ff"  =  0,  /(O)  =  /'(O)  =  0, 

/'(oo)  =  1, 

(2.4) 

and  T]  is  given  by 

1 

(2.5) 
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where  Y*  is  the  Howarth-Dorodnitsyn  variable 

y  =  £  jdY,  (2.6) 

and  we  have  taken  a  model  fluid  (i.e.,  the  Prandtl  number  equal  to  one  and  Chapman’s  vis¬ 
cosity  law  with  C  equal  to  one).  For  the  case  of  a  thermally  insulating  wall,  the  temperature 
is  given  by 

r  =  i  +  5(7-i)Mi(i-/’).  (2.7) 

where  Moo  is  the  Mach  number  of  the  free  stream  and  7  is  the  ratio  of  the  specific  heats  and 
will  be  taken  to  be  1.4  when  it  is  needed  numerically.  Defining  Z  by 
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(?8) 


we  perturb  (2.2)  by  writing 


u  =  Uoo[u  +  8U{X,Y)Eu 


vR-'^  Y8R-'^V{X,Y)Eu  8R-2W{X,Y)Er][l  +  0{R-'^)],  (2.9) 


where  E\  —  exp(iaZ)  and  6  <;  1  (see  Hall  (1982a)  for  further  discussion  of  the  above 
scalings).  We  similarly  perturb  the  basic  pressure,  p,  by  putting 


p  =  p  +  8R-^P{X,Y)Eu 

(2.10) 

and  the  basic  temperature  by 

t  =  Too{T  +  8fiX.,Y)E£, 

(2.11) 

where  Too  is  the  free  stream  temperature.  If  we  then  introduce  a 

,  growth  rate  $  and  the 

scalings 

U{X,Y)  =  U{Y)E2. 

(2.12a) 

V{X,Y)^G'^ViY)E2, 

(2.121)) 

W{X,Y)  =  g'^W{Y)E2, 

(2.12c) 

fiX,Y)^TiY)E2, 

(2.12d) 

P{X,Y)  =  GP{Y)E2, 

(2.12e) 

where 

F2  =exp|y  G'M(A^Kv|, 

(2.12/) 

and  insert  (2.9)  -  (2.11)  into  the  governing  equations  and  take  the  i 

nviscid  limit  G  00  we 

get 

m 

T 

d  fpV\  iap  huT_ 

^dYKT/^T  r  “ 

(2.13a) 

4 


=  -..p, 


rW 

dY' 


u^T  +  1 


,dT 

dY 


0, 


where  K  is  a  measure  of  the  curvature.  Taking  p  =  1  we  find  from  (2.13)  that 


’AV  d  1 

dadu 

d 

hdudY  1 

[t) 

T 

"  W 

du 


+  du 


d  (  1  dV^ 


=  0. 


dY  \TdY^ 

If  we  now  change  variables  to  p,  take  for  simplicity’s  sake  a  value  of  A’  =  1  and  put 


(2.136) 

(2.i3c) 

(2.13d) 

(2.13e) 

(2.14) 


I3=K^3, 


(2.15) 


in  order  to  eliminate  K  from  the  equation  we  get  after  some  rearranging 


V"  +  V 


fff/ 

- - 1- 

/' 


T  f) 


(2.16) 


with  the  boundary  conditions 

1/(0)  =  0,  (2.17a) 

Y  ~  e““''  as  Tj  oo,  (2.176) 

V  ~  — ae~“'’  as  r;  — >  oo.  (2.17c) 

Since  (2.16)  corresponds  to  an  inviscid  limit  we  cannot  satisfy  the  viscous  boundary  condition 
V'  =  0  at  r/  =  0.  Equation  (2.16)  is  the  compressible  generalization  of  equation  (5.8)  in 
Denier,  Hall,  and  Seddougui  (1991).  Similarly  to  Denier,  Hall,  and  Seddougui  (1991)  it  can 
be  seen  that  (2.16)  and  (2.17)  have  the  exact  solution 


=  =  (2.1S) 

which  is  valid  for  all  a  and  Moo-  However,  in  order  to  get  other  solutions  it  is  necessary  to 
solve  this  eigenvalue  problem  numerically. 

We  solved  (2.16)  and  (2.17)  for  a  variety  of  Mach  numbers.  Figure  1  shows  the  first 
eleven  modes  for  a  Mach  number  of  2.  These  modes  are  not  markedly  different  in  shape 
from  those  of  the  form  j3^  =  constant  x  a  obtained  in  the  incompressible  case  (again,  see 
Denier,  Hall,  and  Seddougui  (1991)).  Figure  2  shows  solutions  for  a  Mach  number  of  3.25 
and  we  can  see  that  there  are  now  significant  differences  compared  with  the  incompr^ecibl^* 
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case.  However,  all  of  the  modes  continue  to  rise  as  a  — >■  oo  and  are  of  the  wall  layer  type, 
so-called  because  the  eigenfunction  is  concentrated  in  a  layer  near  the  wall.  In  Figure  3 
where  the  Mach  number  is  now  5,  we  can  clearly  see  the  appearance  of  two  adjustment  layer 
modes  which  tend  to  a  constant  value  of  as  a  — >  oo.  We  call  these  modes  adjustment 
layer  modes  because  their  eigenfunctions  are  concentrated  in  the  temperature  adjustment 
layer  in  the  hypersonic  limit.  (They  will  be  referred  to  as  trapped  layer  modes  for  0(1)  Mach 
numbers.)  Each  of  the  modes  in  this  figure  (apart  from  the  exact  solution)  contributes  a 
portion  to  the  adjustment  layer  modes  as  they  pass  through  the  relevant  positions  and  in 
doing  so  they  change  up  to  the  next  higher  wall  layer  mode.  If  we  follow  the  fourth  mode 
on  Figure  3  as  a  increases  we  can  look  at  its  eigenfunction  and  see  how  it  changes.  Initially 
when  a  —  0.2  we  can  see  in  Figure  4a  that  the  eigenfunction  is  a  widened  version  of  the 
eigenfunction  of  the  third  wall  layer  mode  with  an  additional  lower  peak  in  the  temperature 
adjustment  layer.  In  Figure  4b  with  a  —  1.1  the  mode  is  now  in  its  flat  section  and  the 
eigenfunction  is  indeed  that  of  an  adjustment  layer  mode.  By  the  time  a  =  1.25,  Figure  4c 
shows  that  the  peak  in  the  temperature  adjustment  layer  has  split  into  two.  In  Figure  4d  we 
have  a  =  1.35  and  the  left-hand  peak  of  these  twin  peaks  is  moving  into  the  wall  layer  whilst 
the  right-hand  peak  is  decaying,  showing  the  process  by  which  the  modes  change  up  to  the 
next  higher  wall  layer  mode.  To  complete  the  sequence.  Figure  4e  shows  the  eigenfunction 
for  a  =  2.5  and  it  turns  out  to  be  the  eigenfunction  of  the  fourth  wall  layer  mode. 

Figure  5  shows  solutions  for  a  Mach  number  of  8  and  we  can  see  the  first  three  adjustment 
layer  modes.  As  the  Mach  number  and  wavenumber  increase,  the  modes  come  very  close  to 
intersecting  (although  they  do  not  as  we  shall  show  in  Section  7)  and  this  means  we  have 
to  take  prohibitively  small  steps  in  a  in  order  to  follow  a  particular  mode  numerically.  As  a 
result,  some  wall  layer  type  portions  of  modes  have  been  omitted  from  Figure  5. 


3.  The  Wall  Layer  Modes  for  Moo  '>  1 


Taking  Equation  (2.16)  and  using  the  scalings 

a 


and 


0 


Ml 


J_ 

M« 


-f- 


+ 


(3.1) 


(3.2) 
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allows  us  to  look  at  the  wall  layer  modes  for  Moo  ^  1,  for  w’hich 


J  4/7" 


f2\2  J 


(f‘ 


V"  +  v 


When  a  ^  1  we  can  get  an  approximation  to  the  wall  layer  modes  by  the  application  of  the 
WKB  method.  There  are  three  regions  to  consider  as  we  have  a  turning  point  associated 
with  the  WKB  expansion.  Close  to  the  wall  the  transformation 


V{n)  =  S{im~r\ 


allows  us  to  write  (3.3)  as 


5"(i  -  n  -  25(/7" + /"' + f'n + 5(1  -  n  [(- j«'(7  -  i)'(i  -  ry 

_r  \ ,  K  f'A] 

f  (i-/'2)j^2^^  ;5H/' ^  /)  +  //jj 

which  for  a  :»  1  gives  the  WKB  approximation 

S"  +  a^S  +  =0,  (3.( 


and  we  write  this  as 


S"  +  a^SH{rj,l3)  =  0, 


where 


^  =  -T(7-i)^(i-r)  + 


1(7-1)  (f\,  ,,27  , 


2  B  V  / 


ft  \  ^  ~  J  )  J  J  /■» 


//  >  0  for  0  <  7/  <  T^i  and  r]^  is  the  turning  point.  In  order  tn  satisfy  the  boundary  condition 
(2.17a)  we  must  have  5(0)  =  0  since  (1  —  /')  7^  0  for  r;  =  0.  The  solution  of  (3.7)  which 
satisfies  5(0)  =0  is 

S{v)  =  Asm{af  (//(7/i))2f/r?,),  (3.8) 

H*  Jo 

where  A  is  a  constant  and  so  V  is  given  by 


A{i-n 


i(a  /  (//(7;i))2dr/i). 

Jo 


In  the  region  of  the  turning  point,  the  equation  to  be  solved  is  the  Airy  equation 


V"  -  o^rrU  _  ^  0, 


(3.10) 
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where  a  is  positive  and  given  by 


The  required  solution  of  (3.10)  which  rle<  ays  exponentially  as  ?/  — ♦  cc  is 

I'  =  —  //T). 


(3.12) 


where  B  is  a  constant  and  Ai  denotes  the  Airy  function  (see  Abraniowitz  and  Stegun  (1961)). 
If  we  denote  the  argument  of  the  Airy  function  by  /  then  since  for  /  -C  —  1 


Ai{t)  ^  t:  i  |/|  *  sin  +  — 


we  have  that  for  /  <C  —  1 


n  .  /2 

— r  sin  - 

|/|l  V.3 


Now  letting  tj  — >  //(  in  (3.9)  gives 


A(l-/'(r/,)^)  . 

I  =  - ; - ^Sl 


— ^^sin((i  I  (7/(7/,  ))^(///,). 
—  tl)i  Jo 


77  4  (  7/(  —  7/  )  7  Jo 

Matching  amplitudes  and  phases  between  (3.M)  and  (3.1.3)  implies,  since  we  allow  B  to  be 
positive  or  negative  in  (3.12)  and  hence  phases  may  b('  ('ithi'r  in  phase  or  ~  radians  out  of 
phase,  that 

=±73.  (3.161 


UTT  —  ^ 


where 


./_  =  f  (//(7/i))2  7/;/i 


(3.1S) 


and  71  is  a  large  integer. 

In  Figure  6  we  have  plotted  the  wall  layer  modes,  with  equation  (3.17)  for  n  running 
from  one  through  seven  and  a  Mach  number  of  8,  and  then  su|)('rimposcd  them  on  Figure  5. 
We  can  see,  even  for  a  moderate  Mach  number  and  small  values  of  77,  that  equation  (3.17) 
gives  a  good  approximation  to  the  wall  layer  modes. 


4.  The  Adjustment  Layer  Modes  for  A/^o  1 

For  the  adjustment  layer  modes,  the  appropriate  eqmition  below  the  turning  point  is 
(3.9).  IIow'^‘'''’r.  (3.12),  the  solution  of  the  Airy  equation  (3.10),  must  now  be  renlaced  by 


V  =  Ai{f.)  +  FiTU3i{t), 


(•1.1) 


where  .*7  '  .d  F  are  constants  since  I’  will  be  large  in  the  tenii)erat urt'  acljustinf'nt  layer  away 
from  the  wall  and  so  we  are  no  longer  seeking  a  solution  which  decays  ('Xjjuiumt ially  Ijeyond 
the  turning  point.  Thus  for  ^  <C  —  1,  since 


}3i{t)  ~  '^f\  *  cos  (  j 


we  have 


(£2  ^  f2y,  _  .5.  . 

\  ~  ;  SiU  I  -f-  — ^  —  uj 


lt|t  Vd 

where  tan  o  =  E/F.  Matching  (3.9)  and  (4.3)  gives 

.4(1  -  /'(;;, )^)17<' 


(1.2) 


(4.3) 


± (/•.■*  T  F^F 


(jf' 


and 


where  >i  is  again  a  large  integer. 
When  I  ^  1  we  have 


U7?  -  ^  +  o 


(4.4  i 


(4.01 


and 


so  (4.1)  gives  for  t  3>  1  that 


Bi(t) 


_  1  _  i  /  2  i  \ 

-  -  2/  4  exp  (  -/M  . 


2  s' 

1'  =2  Ft  *  exp  (  -/•'  )  . 


(1.7) 


(4.S) 


and  this  has  to  be  matched  with  the  solution  above  the  turning  ])4)int  but  below  the  tem¬ 
perature  adjustment  layer.  In  this  region  \  is  given  by 

~  f  (_//(, )_  (4.9) 

{  —  If  )*  2^j, 

where  we  have  only  retained  the  exponentially  large  term,  midti|)lied  by  a  constant  cq.  in 
the  \\  KB  approximation.  Matching  betwe«’n  (4.S)  and  (4.9)  gives 


—  =  ^(1  - 
Co  ere 

The  temperature  adjustment  layer  has  c  =0(1)  where 


(4.10) 


rj  =  b+  I 


(4.11) 


and  b  is  a  constant  and 


r  =  /2iog(.\/2 ). 
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In  this  region  therefore 


where  h  is  a  constant  so 


and  hence 


I  ~  1  -f  ('ll  —  l)/(e'. 


(7  -  |(7  -  1)  hen' 


(4.13) 


(4.14) 


(4.15) 


Therefore  if  we  take  2  1  and  insert  these  expansions  into  (4.9)  we  get  that  the  behavior 

r>f  V"  as  the  temperature  adjustment  layer  is  entered  from  below  is 

2co/i2  I  (  a  —  l)lu:' ] 

^  exp]nJ+ - - - — - (4.16) 

Moch-l)^  [  ^  -^‘>0  J 


where 

= r  [i<''  ■ 

In  the  temperature  adjustment  layer  we  make  the  transformation  (4.11)  in  (2.16)  arul 
after  scaling  a  and  l3  by  writing 

a  =  dr  + -  (4.18a) 


0  =  (3r^  +  ..., 

we  get,  retaining  the  leading  order  terms  for  A/oo  1,  that 
cPV  ,  dV  r/l-(7-l)de^\  1  ,  f 


(4.186) 


c^2  ,  1  A5  -  l)/ie‘'\l 


From  (4.16)  we  can  see  that  as  2  ^  00  the  behavior  of  V  will  be  a  multiple  of 


52  e  , 


(4.19) 


(4.20) 


where 


s  =  (7  —  1  )he' . 


Letting  2  — oc  in  (4.19)  we  can  easily  see  that 


1/ 

y  ^  €  , 


since  we  require  V  to  decay  above  the  temperature  adjustment  layer.  Using  (4.20)  and  (4.22) 
as  boundary  conditions  we  solved  (4.19)  numerically.  In  Figure  7  we  have  plotted  the  first 
three  adjustment  layer  modes  given  by  (4.19)  for  a  Mach  number  of  8  and  superimposed 


this  onto  Figure  5.  We  can  see  that  (4.19)  is  giving  a  good  approximation  to  the  adjustment 
layer  modes. 

The  numerical  results  we  gained  in  Section  2  seemed  to  suggest  the  existence  of  a  critical 
Mach  number,  above  which  the  adjustment  layer  modes  come  into  play  and  below  which 
they  do  not  exist,  and  we  shall  discuss  this  further  in  Section  6. 


5.  The  Wall  Layer  Modes  for  Moo  ~  0(1) 


The  modes  located  near  the  wall  for  Moo  ~  0(1)  will  be  described  by  Equation  (2.16). 
For  large  a  we  are  able  to  obtain  an  approximation  to  them  from  the  WKB  solution  in  a 
similar  way  to  the  solution  obtained  in  Section  3  in  the  hypersonic  limit. 

If  we  write 

13  =  I3'^a2 ,  (5.1) 


then  there  is  no  turning  point  in  the  WKB  solution.  However,  the  solution  obtained  does 
breakdown  when  rj  ~  0(<z“*)  as  a  consequence  of  /'  ~  arj  for  7  <C  1.  Thus,  in  this  region  we 
write 

=  aTn,r)^  (5-2) 

where  Tyj  is  the  value  of  T{r])  &t  t]  =  0.  (When  the  flow  is  adiabatic,  i.e.,  the  wall  is  thermally 
insulating,  from  (2.7),  Tu,  =  1  +  (7  —  l)M^/2.)  Substituting  (5.1)  and  (5.2)  into  (2.16)  we 
find  that  V  satisfies  the  following  equation: 


0’ 


(5.3) 


This  is  the  equation  satisfied  by  V  in  the  corresponding  incompressible  problem  described 
by  Denier,  Hall,  and  Seddougui  (1991).  It  is  a  form  of  Whittaker’s  equation  and  the  solution 
for  V  satisfying  the  boundary  conditions 


V  =  0  at  =  0, 


and 


F  — >  0  as  ?/>  — +  00, 


may  be  given  in  terms  of  Kummer  functions  (see  Abramowitz  and  Stegun  1964).  Hence,  we 
find  that  the  unstable  eigenvalues  for  the  present  case  of  a  compressible  fluid  will  be  the  same 
as  those  for  the  incompressible  problem  described  by  Denier,  Hall,  and  Seddougui  (1991). 
These  are 

=  i-L  -1 

’v^’  v/3’-’-  • 
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Thus  from  (5.1)  the  growth  rates  for  the  wall  layer  modes  are  given  by 

=  1^’  (5-4) 

where  n  =  1,2,3,....  We  note  that  the  first  solution  corresponds  to  the  exact  solution 
described  by  (2.18).  Moreover,  the  eigenvalues  are  independent  of  Moo  and  7,„. 

The  appropriate  solution  of  (5.3)  is 

\/(^)  =  -n,2,20),  (5.5) 

where  M{a,b,z)  is  Kummer’s  function  and  6i  is  a  constant.  M{1  —n,2,2rp)  is  a  polynomial 
of  degree  (n  —  1)  in  2'ip.  The  behavior  of  V  for  large  p  is 

K(i/.)~ii(-l)"-‘^e-V  +  -.  -  •  (5.6) 

n! 

Then  from  (2.16),  we  find  that  the  WKB  solution  in  the  region  above  the  wall  layer  which 
matches  with  (5.6)  as  t;  0  is  given  by 

VM  =  c^{T(r|)f-'»HfMrexpl-a  I"  T^dm],  (5.7) 

Jo 

where  cj  is  a  constant.  Since 

/'(?;)  ~  ar;  -  ^7/“  +  ...,  (5.8) 

for  7/  <C  1,  where  a  =  0.332 . . .,  ci  is  given  by 

Cl  =  a"rif"-^>/2a^"(-l)"-^^6i.  (5.9) 

77! 

In  Figure  8  we  show  the  growth  rate  given  by  (5.4)  as  a  function  of  a  for  Moo  ~  5  for 
the  first  nine  wall  modes.  These  solutions  are  superimposed  on  the  corresponding  numerical 
solutions  from  Section  2.  We  see  that  for  large  a  the  asymptotic  results  are  very  good 
approximations  to  the  numerical  solutions. 

In  order  to  discuss  the  near-linking  of  these  wall  layer  modes  and  the  trapped  layer  modes 
described  in  the  next  section,  we  must  consider  the  solution  of  the  wall  layer  modes  for  large 
mode  number.  This  is  apparent  from  Figure  3  since  the  growth  rates  of  the  wall  layer  modes 
will  not  be  close  to  those  of  the  trapped  layer  modes  for  a  1  unless  n  is  large. 

Thus  for  n  >  1  we  write  the  wavenumber  as 

a  =  Gon  +  Oi  +  . . . ,  (5.10) 
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where  Cq  and  oi  are  constants.  Then  from  (5.4)  the  growth  rate  of  the  wall  layer  modes  for 
large  n  is  given  by 


^  =  S  +  ^  +  ---  ■ 

From  (5.5)  the  solution  for  V  in  the  wall  layer  for  large  n  is 

T,,  , .  ,  -_3  _i  _3  i  Stt 

V{'ip)  =  bi2  *Tr  2n  tip*  cos  \8'^ ip — — 


(5.11) 


(5.12) 


As  ^  0  the  solution  (5.12)  matches  onto  the  solution  in  a  thin  layer  next  to  the  wall  of 

thickness  0(n“^)  where  V  satisfies 


SV  2K  ^ 
+  C  “  ’ 


(5.13) 


and  C  =  nil>.  To  satisfy  the  boundary  condition  at  the  wall  and  to  match  with  (5.12)  as 
C  — »•  oo  we  find  that 

V  ~  ei(C  -  +  . . .)  for  C  <  1, 

and 

T,  ,  _i  -1  ..1  /r.1 -i  Stt \ 

V  ~  6i2  *Tr  in  cos  (8242  — —j  as  (—>00, 

where  ei  is  a  constant. 

For  large  values  of  0  the  solution  given  by  (5.12)  is  valid  only  for  ip  <  2n  and  matches 
onto  the  solution  above  the  wall  layer  for  large  n  when  ip  is  large  and  ip  <2n.  We  find  that 
the  WKB  solution  in  the  region  above  the  wall  layer  which  matches  with  (5.12)  as  77  0  is 

given  by 

y{v)  =  gi{-Q{T]))~iT{r])  cos{{nao  +  ai)  j  {-Q)^dr]i-^  f  I  ■  ^  +  (-Q)A  dgi-9i), 

Jo  Zao  Jo  \(_(5)2  J 

(5.14) 

where  gi  is  a  constant  and  61  =  Ztt/A.  The  function  Q{g)  is  defined  by 


Qiv)  =  (n^)r  - 


,  2  (ng)f"{g)  T{g) 


(5.15) 


and  the  solution  (5.14)  is  valid  for  0  <  77  <  77*  where  77*  is  the  first  zero  of  Q  and  Q  <  0  for 
77  <  77*.  A  more  thorough  discussion  of  the  function  Q  is  given  in  the  next  section. 

In  the  region  close  to  77  =  77*,  from  (2.16)  V  satisfies 


V"  -  T{g  -  g)aln:^V  = 


(5.16) 


where  r  =  Q'{g‘)  >  0.  Since  we  require  V  to  decay  as  77  —7  00  the  appropriate  solution  of 
(5.16)  is 

V  =  koT^iAi{r),  (5.17) 
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where  ko  is  a  constant  and  r  =  (con^r)^  (77  —  77*).  From  the  behavior  of  the  Airy  function  Ai 
given  by  (3.13)  for  r  —00  we  have  from  (5.17)  that 


V  ~  kolrl  <  cos  > 

(5.18) 

as  r  — >  —00.  Then  matching  the  amplitudes  and  phases  of  (5.14)  as  7/ 

— +  Tj*  with  (5.18)  gives 

'r~Haln‘^T)^T{j]*)gi  =  ko, 

(5.19a) 

and 

(nao  +  ai)Ii  —  - — (Ii  -f  I2)  —  0i  = 

Zqq  4 

i5.19b) 

where  the  integrals  7i  and  I2  are  defined  by 

h  =  1  {-Q)^dT)i, 

Jo 

(5.20a) 

and 

h=  T\-Q)-Ur),. 

Jo 

(5.206) 

The  expression  (5.19b)  determines  aj,  i.e., 

+  4  ~  naoli 

(5.21) 

h-^,{h+hy 

6.  The  Trapped  Layer  Mode  for  Mo 


For  Moo  ~  0(1)  this  mode  is  concentrated  away  from  the  wall  and  is  the  equivalent  of 
the  temperature  adjustment  layer  made  for  Moo  ^  1  described  in  Section  4.  For  a  >■  1  this 
mode  can  be  regarded  as  a  virtually  continuous  function  of  a.  In  Section  7.2,  we  discuss  the 
near-linking  of  this  mode  with  the  wall  layer  solution  for  large  n.  Anticipating  this  we  define 
a  by  (5.10)  and  consider  the  solution  of  the  trapped  layer  mode  for  n  »  1.  To  obtain  the 
solution  for  large  a  from  the  following  results  we  set  oi  =  0  and  non  =  a.  For  large  n  we 
write  (2.16)  in  the  form 


where 


V"  -2^V'  -aln^QV  = 


Then  this  mode  is  concentrated  in  the  region  where 


= 


Tf  2T 
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and  /3  is  a  maximum  bX  rj  =  fj.  Therefore,  from  (6.2)  we  have 

Q{.V)  =  Q'iv)  =  0  and  Q"{n)  >  0.  (6.4) 

Thus  when  cq  =  2/?o  the  function  Q  defined  by  (6.2)  is  also  that  defined  by  (5.15). 

In  Figure  9  we  plot  Q{ri)fT  from  (6.2)  for  A/qq  =  5  for  a  range  of  values  of  (3.  We  choose 
to  show  Q IT^  ather  than  Q  since  {T{t])Y  given  by  (2.7)  is  large  for  small  values  of  rj  and 
decays  exponentially  to  the  free  stream  value  of  1  as  7?  increases.  We  see  from  Figure  9  that 
Qi'n)  BlS  T]  0+.  For  certain  values  of  ^  the  function  Qir])  has  three  zeros.  For  a 

particular  value  of  ^  the  second  and  third  zeros  coincide  in  a  turning  point.  This  position  of 
■q  is  the  location  of  the  trapped  layer  modes  q  =  q,  for  large  a.  The  sequence  of  modes  will 
be  described  by  corrections  to  (6.3). 

However,  these  trapped  modes  will  not  exist  for  values  of  below  a  critical  value.  We 
can  see  this  if  we  look  at  the  plots  of  Q{q)jT  for  various  values  of  (3  for  M^o  =  3  shown  in 
Figure  10.  It  is  apparent  from  this  figure  that  Q{q)  only  has  one  zero  for  any  value  of  (3  and 
no  local  minimum  for  A/oo  =  3.  For  the  Chapman  constant  C  =  1  and  an  adiabatic  fluid  the 
critical  value  of  Moo  is  found  to  be  =  Me  —  3.564.  Thus,  for  Moo  <  Me  the  function  Q 
will  not  have  a  local  minimum  and  there  is  no  solution  described  by  (6.3)  and  (6.4).  This 
explains  why  the  discontinuous  modes  do  not  appear  in  the  computational  results  for  small 
Mach  numbers.  We  note  that  there  exists  a  critical  value  of  the  free  stream  Mach  number 
when  C  7^  1  or  Prandtl  number  ^  1  and  also  for  the  case  of  an  isothermal  flow  where  the 
temperature  of  the  wall  is  maintained  at  a  constant  value. 

Thus,  from  (6.3)  /?  is  a  constant  to  first  order  for  the  trapped  modes  and  moreover,  the 
same  constant  for  each  mode.  The  correction  to  ^  will  describe  the  distinct  modes.  In  order 
to  determine  the  eigenfunction  in  this  region  it  is  necessary  to  determine  the  correction  to 
(3.  In  the  trapped  layer 

i  1 

q  =  a^n^iq  -  q),  (6.5) 

where  q  0(1).  Write 

/3  =  Po  +  cio^n  ^^1  +  ...,  (6.6) 

and  substitute  (6.5)  and  (6.6)  into  (6.1).  Equating  coefficients  of  the  largest  terms,  which 
are  of  O(aQn^),  shows  that  /?o  is  indeed  defined  by  (6.3)  and  (6.4). 

Figures  11  and  12  show  the  solutions  of  (6.3)  and  (6.4)  for  an  adiabatic  fluid  and  also  for 
the  case  of  an  isothermal  fluid  where  the  temperature  of  the  wall  is  maintained  at  a  constant 
value.  In  this  case  for  Chapman  constant  C  =  1  and  a  Prandtl  number  of  unity,  the  basic 
temperature  is  given  by 

nv)  =  1  +  —MKnri)  -  /'(,))>)  +  (r.  -  1)(1  -  /'(,)),  (6.7) 


15 


where  Tyj  is  the  nondimensional  temperature  of  the  wall.  Note  that  if  Tu;  =  1  +  (7  —  l)M^/2 
then  (6.7)  reduces  to  the  expression  (2.7)  for  an  adiabatic  flow. 

Figure  11  shows  the  value  of  rj  for  the  trapped  modes  as  a  function  of  Mach  number 
for  an  adiabatic  wall  and  an  isothermal  wall  with  =  0.5,  1,  2,  4,  8.  VVe  see  that  fj 
moves  towards  the  free  stream  as  Moo  increases.  For  large  values  of  Moo  the  analysis  must 
be  replaced  with  that  described  in  Section  4  so  7  will  correspond  to  the  position  of  the 
temperature  adjustment  layer  given  by  (4.11)  and  the  modes  will  be  as  described  in  Section 
4.  This  is  consistent  with  (4.11)  since  we  see  that  the  position  of  the  temperature  adjustment 
layer  increases  as  Moo  increases. 

Figure  12  shows  the  corresponding  values  of  /3o  for  the  trapped  modes  as  a  function  of 
Mach  number.  We  see  for  the  adiabatic  case  and  for  the  isothermal  case  with  <  S  that  the 
value  of  /3o  decreases  initially  as  Moo  increases  but  then  increases  as  Moo  increases  further. 

The  O(aon)  terms  in  (6.1)  give  the  following  equation  for  V: 


(PV 

drf' 


Q"{fi) 

2 


ev  -  ‘^(T(fi)fv  =  0. 


(6.8) 


where  Q"{rj)  is  defined  for  f3  =  /?o.  The  solution  of  (6.8)  satisfying  the  boundary  conditions 


V  -+  0  as  I?)!  — >  00, 


(6.9) 


is  an  eigenvalue  problem  for  We  find  that  the  solution  of  (6.8)  and  (6.9)  is 


where  di  is  a  constant, 


Virj)  =  drU{X,fi), 

(6.10) 

V  =  (2Q"(»)))H, 

(6.11) 

^  22(r(7))2^i 

(Q"(77))2/?o  ’ 

(6.12) 

and  U{X,T))  is  the  Parabolic  cylinder  function  (see  Abramowitz  and  Stegun  1964).  Now 
U{X,r))  grows  exponentially  as  rj  —00  unless 


A  = 


1  _3  5 

2’  ”2’  ~2’‘"  ■ 


(6.13) 


Then  the  behavior  of  U{X,t])  as  7/  — >  —00  when  A  is  given  by  (6.13)  is 


and 


f/(A,77)~(-l)'^-^l7|-^-^e-^^^ 

(6.14) 

(6.15) 
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as  7]  +00. 

From  (6.12)  and  (6.13)  we  can  evaluate  /3i  for  a  fixed  iVIach  number  for  a  particular 
trapped  mode.  Figure  13  shows  the  growth  rate  from  (6.6)  for  Moo  —  5  for  the  first  three 
trapped  modes  as  a  function  of  a.  In  order  to  see  the  comparison  between  the  numerical 
solution  of  the  discontinuous  mode  described  in  Section  2,  in  Figure  13  only  the  discontinuous 
solutions  are  shown  for  large  a  and  below  a  certain  value.  We  see  from  Figure  13  that  for 
the  first  trapped  layer  mode  the  asymptotic  solution  is  a  very  good  agreement  even  for  0(1) 
values  of  a.  For  successive  modes  we  still  have  good  agreement  but  for  larger  values  of  a. 

Below  T]  —  f}^  but  above  r]  =  t}*  where  Q{r)‘)  =  0,  Q{ri)  will  be  positive.  The  WKB 
solution  of  (2.16)  in  the  region  rj’  <  rj  <  fj  is 


y{v)  =  ^o{Q{r]))  *T{7j)exp 


'  4. 

aoTi  +  ai  -  — 
i  Po 


(6.16) 


where  cq  is  a  constant.  Now  as  r/  — > 


(6.17) 


and 


jSi  r>7-(aon)“^ 


T'Q-5d7;i  -  A  ln{{2Q"{fi)r^rj\),  (6.18) 


where  A  satisfies  (6.12).  Thus,  from  (6.16)  as  r/  — ^  for  large  n 


f 

Jtj* 


V~2ta^ni(Q"(77))-«r(7))eoM-"-^(2Q")"exp 


+  (  +  Cl 


I) 


,  (6.19) 


where  I3  and  are  integrals  defined  by 


h  = 


(6.20a) 


and 


Pa  Jti* 


Then  (6.19)  matches  with  (6.10)  as  r/  -00  if 


2hlnHQ"m-^T{f,){2Q")^eoexp 


01 


aoTi  +  Oi  —  —  /a  +  /4 

V  PoJ 


=  di(-l) 


(6.206) 


(6.21) 


Close  to  the  position  rj  =  t}*  V  satisfies  (5.16).  We  do  not  require  V  to  decay  above  this 
region  so  the  solution  for  V  is 


V  =  koTT^  Ai{r) Iqtt^  Bi{r), 


(6.22) 
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where  ko  and  Iq  are  constants. 

From  (4.8)  the  behavior  of  K  as  r  — >  oo  is 


—  —  -N 

V  lor  *  exp  (  )  • 


The  solution  (6.16)  must  match  with  (6.23)  a.s  j]  rj*.  Then  /q  is  given  by 


lo  =  eoT  4(a^n^r)i2T(7/*). 


From  (6.22)  as  r  — >  — oo 


^  ~  (^0  +  *  cos  ^  -  02^  , 


(6.23) 


(6.24) 


(S.zdj 


where  tan  02  =  ko/lo.  In  the  region  below  Tf  =  tj*  V  is  given  by  (5.14),  i.e.,  the  solution  above 
the  wall  layer  for  n  >■  1.  Thus  matching  the  solution  (5.14)  as  t;  — >  77*  with  (6.25)  gives 
instead  of  (5.19) 

T-*(a^n^T)^T(Tf*)gi  =  ±(k^  +  l^)^,  (6.26) 


(nao  +  ai)Ii  —  - — (Ii  +  I2)  —  ^1  —  ^2  “  ■j- 
Zcto  4 


(6.27) 


7.  The  Near  Crossing  of  the  Wall  Layer  and  Adjustment  Layer  Modes 

We  see  from  the  numerical  results  of  Section  2  that  for  values  of  the  free  stream  Mach 
number  above  a  critical  0(1)  value  the  wall  layer  modes  and  the  adjustment  layer  modes 
become  very  close.  This  occurs  for  0(1)  values  of  the  wavenumber  and  continues  to  occur  eis 
the  wavenumber  increases. 

A  similar  near  crossing  occurs  for  the  inflectional  acoustic  neutral  modes  and  the  vorticity 
mode  associated  with  the  inviscid  instability  of  a  Blasius  boundary  layer  for  a  compressible 
fluid  and  was  discussed  by  Cowley  and  Hall  (1990)  and  Smith  and  Brown  (1990).  The  latter 
authors  confirmed  the  conjecture  of  the  former  that,  in  the  hypersonic  limit,  the  modes  are 
separated  by  an  exponentially  small  amount. 

To  discuss  the  near  crossing  of  the  two  different  types  of  mode  described  in  Sections  3 
and  4  for  Moo  ^  1  ^^nd  in  Sections  5  and  6  for  a  1  we  follow  Smith  and  Brown  (1990) 
and  extend  our  WKB  analysis  of  these  previous  sections. 

Since  the  near  crossing  occurs  for  0(1)  Mach  numbers  as  well  as  for  Moo  ^  1  it  should 
be  possible  to  discuss  the  near  crossing  for  both  situations.  The  analvsis  will  be  similar  for 
each  case  and  therefore  we  attempt  to  minimize  any  repetition.  We  will  first  discuss  the 
situation  for  the  hypersonic  limit. 


7.1.  The  Hypersonic  Limit 


We  again  have  that  below  the  turning  point  V  is  given  by  (3.9).  In  the  region  of  the 
turning  point  V  is  given  by  (4.1)  and  (4.3)  -  (4.5)  still  hold.  However,  we  replace  (4.8)  by 

V  fs  ^  ^jEexp  +  Fexp  ,  (7.1) 

for  <  ^  1.  We  also  retain  the  exponentially  small  term  in  (4.9)  so  that  V  in  the  region 
between  the  turning  point  and  the  temperature  adjustment  layer  is  given  by 

coexpla/  (-H(T]i))^dT]i\  +  doexp  l-a  f  {  -  H  drjA  .  (7.2) 

( —  H)*  t  Jilt  )  L  Jrit  > -i 

We  already  have  (4.10)  and  matching  the  exponentially  small  terms  between  (7.1)  and  (7.2) 


gives 


=  ^(1  - 


4 = -■  (7-^) 

,  2do  cq 

We  now  need  to  augment  (4.16),  the  behavior  as  the  temperature  adjustment  layer  is  ap¬ 
proached  from  below,  by  a  corresponding  term  in  do  so  that 


M<«(7  -  1) 


^  1,  r  f  _  T  a  (7  —  l)h  ,  f  _T  a  (7  —  l)h  , 

-^62"  |^Coexp|aJ+  -  -  |  do  exp  |-a7+  -h  -  e 


and  then  match  this  with  the  behavior  of  V  in  the  temperature  adjustment  layer  as  2  — +  00. 
If  we  assume  (4.19)  has  a  solution 


a  =  A, 
0  = 


(7.6a) 

(7.66) 

(7.6c) 


and  then  perturb  these  equations  by  putting 


a  —  /Id"  v4. 


(7.7a) 

(7.76) 


P=B  +  B, 


(7.7c) 
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and  insert.  (7.7)  into  (4.19)  we  get  after  linearising 


=  V  2AA(1 -Kj  -  l)he‘y  + 


A^B  AA 


(7  —  \)he^ 

2 


which  can  be  integrated  to  give 


V  =  Vw{z), 


w{z)  =  2AA 


/.'( 


,-r3  r /[i  +  (7-l)/ie^'P  r'T72 


y  d22  dzi 


^  ([^  +  il  -  l)/ie''*]^ 


Cl  n  [1  +  (7  -  l)fee^>]^ 
(7  -  l)h  Jo 


[1  +  (7  - 


dzi  -|-  C2. 


dz2  dzi  (7.10) 


Using  (4.22)  we  can  calculate  the  constant  Ci  that  is  necessary  to  ensure  that  V  decays 
exponentially  as  z  — ♦  —00.  We  can  then  use  (4.20)  to  calculate  the  behavior  of  w{z)  as 
z  —*  00  and  if  we  choose  C2  so  that  there  is  no  additional  contribution  of  order  S2e““*  to  V 
from  V  then  we  have  that  as  2  — >  00,  U  goes  as 


(7.11) 


^  ~  [(7  ~  l)he^Y  k-'^p{~-^(7  “  l)^e^}  + 


V  B^ 


AB  A  \  —  l)h 


Ca  exp{y4(7  —  l)/ie*}  , 


where 


V^dz 


V'^dz, 


(7.12) 


roo 


^  Jo  [1^  +  (7  “  l)/ie*]^  Jo  ] 

Matching  between  (7.5)  and  (7.11)  gives  us  that 

^  =  -Uc3+fl|-4)^^ 

do  \B  B  2 


[1  +  (7  -  \)he^Y 


Ca  exp{— 2aJ+}, 


(7.13) 


(7.14) 


20 


so  we  can  get  Fj E  from  (7.4).  Assuming  F  is  small  since  we  want  to  perturb  the  wall  layer 
modes  (as  well  as  the  adjustment  layer  modes)  for  which  F  is  zero,  we  get 


(7.15) 


and  inserting  this  into  (4.5)  gives  after  some  rearranging 


mr  —  j 

Mf- 


AB{i -J)hC4  ]\ 
{-f  -l)hC4B  -2C3B^\) 


_ B^T _  exp{-2aM^J+} 

2C3B^  -  (7  -  l)hC4  MIJ. 


(7.16) 


Setting  the  two  brackets  equal  to  zero  gives  respectively  the  wall  layer  modes  and  the  adjust¬ 
ment  layer  modes.  We  can  see  from  (7.16)  that  for  finite  Moo  these  modes  do  not  intersect 
but  the  distance  between  them  is  exponentially  small.  This  explains  the  difficulties  we  had 
in  following  a  particular  mode  when  we  solved  (2.16)  -  (2.17)  numerically.  A  sketch  is  given 
in  Figure  14  of  a  pair  of  these  near  crossings  with  the  continuous  lines  representing  the  situ¬ 
ation  at  a  finite  Mach  number.  The  dashed  lines  represent  the  picture  given  by  (7.16)  with 
the  exponentially  small  right-hand  side  ignored,  i.e.,  when  the  adjustment  layer  modes  have 
become  continuous  in  the  limit  Moo  oo. 


7.2.  Moo  0(1) 

As  discussed  in  Section  5,  in  order  to  describe  the  near-linking  of  the  wall  layer  modes 
and  the  trapped  layer  modes  for  Moo  ~  0(1)  we  must  consider  the  situation  when  the  mode 
number  of  the  wall  layer  mode  is  large.  Similarly  to  the  case  when  M^o  ^  1  we  extend  the 
WKB  analysis  of  Sections  5  and  6  for  n  >>  1. 

The  solution  in  the  wall  layer  for  V  is  given  by  (5.12)  with  the  solution  below  the  position 
T]  =  rj*  given  by  (5.14).  The  solution  close  io  t}  —  r]*  is  given  by  (6.22)  and  so  (6.26)  and 
(6.27)  are  still  satisfied. 

Now  as  r  — >  oo,  from  (6.22)  we  replace  (6.23)  by 

_i  ko  /  2  3\  /2  sN 

Vr^r  * 

In  the  region  tj*  <  r]  <  fj  we  must  retain  the  exponentially  decaying  term  so  instead  of  (6.16) 
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we  have 


y{v)  =  (Qiv))  ‘T(7/)|eoexp  ^nao  +  ^  +  ~ 

+  /(oexp  -  (^nao  +  cii  -  J ^  Q^drji  “  ^  |  , 

(7.18) 

where  ho  is  a  constanl.  Matching  this  solution  as  rj  ^  rj*  with  (7.17)  gives  (6.24)  as  well  as 


A’o  _  ho 

2/o  Co 

Following  (6.19),  from  (7.18)  as  rj  —y  fj  for  n  large 


(7.19) 


V  ~  28a^nf(C?"(7}))  *'^(^)  {  eo|//|  ^  2(2(5"(j7))<"  exp  -'-r  +  [  +  ai  -  ^  h  +  h 


+/iol7/l'^  '^{2Q''{fi))  "'exp  ~  -  [aon  +  ax  -  I3  -  I4  L 


(7.20) 

Now  we  make  the  transformation  (6.0)  and  perturb  the  solution  in  the  trapped  layer  lor 
large  n  where  a  is  given  by  (5.10).  Thus  we  write 


^  —  /^o  4 - (/?!  +  ^*)  +  .  • . : 

aon 


(7.21) 


V  =  Fo  +  K*, 


(7.22) 


where  /?*  and  V*  are  small  perturbations  and  Vo  satisfies  (6.8).  Then  if  we  substitute  (7.21) 
and  (7.22)  into  (2.16)  and  linearize  about  the  perturbed  quantities  we  find  that  V*  satisfies 


(7.23) 


We  wish  to  determine  the  behavior  of  K  =  Vq  +  V*  as  77  — >  — cx3.  We  want  ^‘(t))  to 
decay  as  77  —>  00  so  that  the  disturbance  does  not  propagate  outside  the  boundary  layer. 
The  solution  of  (7.23)  is 

V(^)  =  ♦(i))V'o(^),  (7.24) 

where  is  given  by 


’’  drjx 
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Using  the  behavior  of  Vo  as  77  ±00  given  by  (6.15)  and  (6.14)  we  can  determine  the  values 

of  the  constants  Di  and  D2-  Then  from  (6.15),  since  V'  must  decay  for  large  7}  we  find 

9/1*  _  1  /•<» 

D,  =  {UiXM'^dJ}.  (7.26) 

For  T/  — >  —00  we  choose  D2  to  cancel  the  constant  terms.  Thus,  from  (6.14)  we  find  that 

^(W)  ~  T>3l^P'exp(7/^/2),  (7.27) 


as  7/  — >  —  00  where 


T>3  = 

1  F'?/!’  r-oo  , 

/  (t/(A,77))^d^(2Q"(77))-  + 

Po  2o 

A' 

di 

(7.28) 

Thus  from  (7, 

.26), (7.28) 

is 

4/i*  , 

Ds  = 

Po 

(7.29a) 

where 

h=  r(U(K'i)fdJi- 

Jo 

(7.296) 

Thus,  from  (1 

’.22),  (7.24) 

,  (7.27),  and  (6.14),  the  behavior  of  V  as  ij  — > 

—  (OC' 

is  given  by 

1^(77)  ~  (fi(-l)"'^  2  e.xp(-r7V4)  +  'Pxp(77V4)].  (7.30) 


The  VVKB  solution  below  the  trapped  layer  must  match  with  (7.30)  as  77  — >  —00.  Thus 
matching  (7.20)  and  (7.30)  gives 


Co  {2Q"{fi))  2  /I, 

=  - =; - exp[-2(nao  +  a,  -  — )/;5  -  27,,]. 

Ito  L>3  yo 


(7.31) 


Thus,  from  (7.19)  we  have  /q/^'o-  The  .solution  for  V  described  above  is  also  a  perturbation 
of  the  wall  layer  solution  so  we  must  have  from  (6.22)  that  /q  is  small.  Then  we  have 


TT  /o 

2  ko 


(7.32) 


Substituting  this  in  (6.27)  gives 


(h{h  —  ~ — {h  +  h))  -k  naoI\  —  0\  -  —  — — . 

2ao  4  K\) 


(7.33) 


Using  (7.19)  and  (7.31)  we  can  rearrange  (7.33)  to  give 

O1  +  ^  -  naoU  W{l))-2  exn[-2(7iuo  +  e,  -  -  2/.,] 

■  A  -  i(7i  +  7«)  277,(7,  -  +  /,)) 


(7.3-1) 
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Now  from  (5.11) 


ax 


8  2  Aq  ^2 


(7.35) 


Substituting  the  expression  for  gi  from  (7.34)  into  (7.35)  and  rearranging  gives 


1  /0i  4-  7r/4  —  noo/i 


SsGqG 


.  +  ^2) 


0  -  00  - 


A 

uon 


(7.36) 


Sla!7r{T{^)yh{h-Mh  +  h)) 


exp 


0x 


-2  I  nao  +  ^  1  /s  -  21. 


If  we  set  the  factors  on  the  left-hand  side  of  (7.36)  to  zero  in  turn  we  see  that  they 
describe  the  growth  rates  of  the  wall  layer  modes  and  the  trapped  layer  modes  respectively. 
Thus  from  (7.36)  we  see  that  for  large  n  the  wall  layer  modes  and  the  trapped  layer  modes 
are  separated  by  an  exponentially  small  amount.  This  explains  the  apparent  crossings  of 
the  growth  rate  curves  obtained  fmm  the  numerical  results  of  Section  2.  Thus  the  curve  of 
growth  rate  as  a  function  of  a  where  the  modes  are  very  close  is  also  described  by  Figure  14 
as  for  the  case  when  Moo  ^  1- 


8.  Discussion 

The  main  result  of  the  present  paper  is  that  we  are  able  to  investigate  the  linear  growth 
rate  of  Gortler  vortices  for  0(1)  wavenumbers.  This  was  achieved  by  considering  the  inviscid 
limit  of  large  Gortler  number  for  vortices  having  O(G^)  growth  rates.  In  this  limit  the 
growth  of  Gortler  vortices  is  governed  by  parallel  flow  effects.  Since  for  incompressible  flows, 
or  compressible  flows  with  0(1)  Mach  numbers,  the  spatial  growth  of  viscous  Gortler  vortices 
is  governed  by  non-parallel  effects  previous  investigators  were  forced  to  consider  the  limit 
of  large  wavenumber,  where  as  a  result  of  boundary  layer  growth,  non-parallel  effects  are 
unimportant  (sec  Hall  1982a  and  Hall  and  Malik  1987).  However,  Hall  and  Fu  (1989)  show 
that  in  the  hypersonic  limit  the  growth  of  Gortler  vortices  with  wavelength  0((2Iog  4/^  )“? ) 
is  governed  by  a  parallel  flow  theory. 

We  have  shown  that  in  the  inviscid  limit,  in  a  compressible  boundary  layer  over  an  infinite 
cylinder,  there  are  two  types  of  growth  rate  modes  possible  for  Gortler  vortices. 

In  the  hypersonic  limit  we  have  firstly  the  type  of  mode  we  call  wall  layer  modes.  These 
are  [)resent  m  the  incompressible  case  and  their  growth  rates  continue  to  rise  as  the  wavenum- 
b(>r  a  — +  00.  Secondly,  we  have  the  adjustment  layer  modes  which  have  their  eigenfunctions 
concentrated  in  the  temperature  adjustment  layer  away  from  the  wall.  The  growth  rate  of 
these  modes  tends  to  a  constant  value  as  a  — >  00. 
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We  have  also  considered  the  limit  a  — >  oo  for  0(1)  Mach  numbers  to  describe  the  wall 
layer  modes  and  the  trapped  layer  modes  (equivalent  to  the  temperature  adjustment  layer 
modes  when  M^o  !)•  The  wall  layer  modes  are  confined  to  a  thin  layer  of  thickness  0(a~* ) 
and  have  growth  rates  proportional  to  02  for  moderate  mode  number  n  and  tending  to  a 
constant  value  for  large  n.  These  results  are  the  extension  for  compressible  flows  of  the 
incompressible  case  described  in  Section  5  of  Denier,  Hall,  and  Seddougui  (1991),  hereafter 
referred  to  as  DIIS.  For  the  compressible  case  we  have  an  additional  mode  concentrated  in 
a  layer  of  thickness  0(a~2)  away  from  the  wall.  This  is  the  so-called  ‘‘trapped  layer  mode,” 
which  has  growth  rates  tending  to  a  constant  value  as  a  — >  00. 

The  asymptotic  results  described  for  Moo  1  and  also  for  Moo  ~  0(1)  are  shown  to 
agree  very  well  with  the  numerical  solutions  described  in  Section  2.  We  find  that  as  the 
Mach  number  is  increased  from  zero  the  wall  layer  modes  start  to  deform  until  at  a  critical 
Mach  number  we  see  the  appearance  of  adjustment  layer  modes. 

These  modes  are  not  present  in  an  incompressible  fluid.  However,  the  growth  rate  of  these 
mo-  .s  tends  to  a  constant  value  as  the  wavenumber  increases  while  that  of  the  wall  layer 
modes  (also  present  in  incompressible  flows)  continues  to  increase.  Hence,  we  anticipate  that 
the  trapped  layer  modes  will  not  be  as  important  as  the  wall  layer  modes  for  large  values  of 
the  wavenumber. 

The  situation  here  is  similar  to  the  case  of  inviscid  disturbances  to  coni pi'essible  flow  over 
a  flat  plate.  In  this  instance  the  vorticity  mode  does  not  appear  until  Moo  exceeds  a  value 
approximately  given  by  2.2.  However,  in  contrast  to  the  present  problem,  the  vorticity  mode 
has  larger  growth  rates  than  the  acoustic  modes  when  A/co  ^  1  • 

From  the  analysis  for  a  1  for  Chapman  constant  C  =  1  and  unit  Prandtl  number 
for  an  adiabatic  boundary  condition  on  the  basic  temperature  we  find  that  this  critical 
Mach  number  is  3.564.  For  Mach  numbers  of  this  range,  the  adjustment  layer  modes  are 
discontinuous  with  each  mode  solution  of  (2.16),  apart  from  the  exact  solution,  contributing 
a  part  to  the  adjustment  layer  modes  as  it  passes  through  the  relevant  positions. 

We  showed  in  Section  7.1  that  for  Moo  ^  1  the  wall  layer  solutions  and  temperature 
adjustment  layer  solutions  of  (2.16)  come  v  ithin  an  exponentially  small  distance  of  one 
another  and  in  the  hypersonic  limit  the  adjustment  layer  modes  become  continuous.  In 
Section  7.2  we  showed  that  this  is  also  the  situation  in  the  limit  a  —>  00  for  Moo  ~  0(1)  when 
the  mode  number  of  the  wall  layer  modes  is  large. 

Idle  Chapman  viscosity  law  has  been  assumed  in  the  present  analysis,  and  results  dc- 
scril>ed  for  unit  Chapman  constant  and  unit  Prandtl  number  for  simplicity.  There  is  no 
great  diflicidty  in  obtaining  results  for  more  realistic  values  of  C  and  Prandtl  number  and  it 
is  expected  tliat  tlu'  effect  of  this  on  the  asymptotics  will  be  small. 
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However,  in  the  hypersonic  limit  the  Chapman  viscosity  law  does  not  gi’ a  realistic 
description  of  the  viscosity  of  the  fluid.  Thus,  of  interest  would  be  the  results  when  the 
more  accurate  relation  Sutherland’s  law  is  used  to  describe  the  temperature  dependence 
of  the  fluid  viscosity.  The  results  of  Fu,  Hall,  and  Blackaby  (1990)  who  investigated  the 
Gortler  instability  in  the  hypersonic  limit  for  Sutherland’s  law  showed  that  the  results  were 
significantly  different  from  those  for  Chapman’s  law.  In  particular  the  structure  of  the 
temperature  adjustment  layer  is  altered,  now  having  0(1)  thickness.  See  also,  Blackaby, 
Cowley,  and  Hall  (1990)  who  investigated  hypersonic  flows  over  a  flat  plate  using  Sutherland’s 
law.  Thus,  it  is  to  be  expected  that  significant  changes  to  the  results  presented  in  this  paper 
would  exist  if  Sutherland’s  viscosity  law  was  used  in  place  of  Chapman’s  law. 

The  solutions  of  the  equations  governing  the  instability  of  Gortler  vortices  in  the  inviscid 
limit  do  not  predict  a  fastest  growing  mode  since  the  growth  rates  of  the  wall  layer  modes 
tend  to  infinity  when  a  oo.  From  Section  5,  in  the  high  wavenumber  limit  we  see  that  the 
wall  layer  modes  have  growth  rates  of  0(G2<Z2).  We  can  extend  the  work  of  DHS  to  describe 
the  fastest  growing  mode  in  the  limit  of  large  Gortler  number  for  a  compressible  flow  with 
Moo  ~  0(1).  First  we  consider  the  structure  of  a  viscous  mode  close  to  the  right-hand  branch 
of  the  neutral  curve  with  a  =  X*G*.  This  mode  exists  in  a  layer  of  thickness  0(a“5)  centered 
around  some  non-zero  value  of  t].  Initially  A*  is  0(1)  but  we  are  interested  in  the  limit  A*  — +  0. 
We  note  that  the  problem  where  A*  — >  oo  corresponds  to  the  work  of  Hall  and  Malik  (1989). 
The  analysis  is  very  similar  to  that  described  by  DHS  for  the  corresponding  problem  for  an 
incompressible  problem  so  the  details  will  not  be  repeated  here.  As  for  the  inviscid  modes 
the  growth  rates  are  0(G2).  The  first  approximation  to  the  growth  rate  4  is  determined 
from  the  following  condition 


Ku  uy 
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=  0,  (8.1) 


and  ^  is  a  maximum  at  y  =  yc  where  (8.1)  is  satisfied.  Here  p  is  the  viscosity  of  the 
fluid  and  F*  is  the  Prandtl  number.  In  the  incompressible  limit  (8.1)  reduces  to  precisely 
the  expression  given  by  DHS  for  the  corresponding  incompressible  problem.  We  consider 
solutions  of  (8.1)  when  the  basic  flow  is  given  by  the  compressible  Blasius  flow  described  by 
(2.4  -  7)  for  unit  Chapman  constant  and  Prandtl  number.  Then  the  situation  of  A*  =  0  is 
precisely  (6.3),  determining  the  growth  rate  of  the  trapped  layer  mode.  Thus,  for  A*  =  0 
there  will  not  be  a  solution  of  (8.1)  for  values  of  Moo  below  a  certain  value.  However,  for 
A*  >  0  we  find  that  solutions  of  (8.’)  exist  for  all  values  of  Moo-  We  are  not  concerned  with 
the  details  of  the  solution  of  (8.1)  here  but  further  particulars  may  be  obtained,  on  request, 
from  the  second  author.  We  are  solely  concerned  with  the  behavior  of  the  solution  of  (8.1) 
for  A*  0  for  the  basic  flow  described  in  Section  2.  The  numerical  solutions  of  (8.1)  show 
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that  T)c  becomes  small  for  A*  — >  0  and  closer  inspection  of  (8.1)  in  this  limit  reveals  that 
r]c  ~  A*‘*  and  /3  ~  A*“^.  Thus,  in  the  limit  A*  — )•  0  the  growth  rates  of  the  viscous  modes  are 
0(G2(a(j~<  )~^).  Thus,  as  a  decreases  the  viscous  right-hand  branch  modes  become  more 
unstable. 

The  above  discussion  suggests  that  there  exists  an  intermediate  region  where  the  viscous 
mode  described  above  and  the  inviscid  mode  of  Section  6  overlap.  This  will  occur  when 
their  growth  rates  are  the  same  size.  We  find  that  this  is  the  case  when  a  ~  with  the 

3 

growth  rates  of  0(G5).  The  vortices  in  this  case  will  be  confined  to  a  wall  layer  of  depth  of 
0(G'”s).  We  note  that  this  situation  is  identical  to  that  described  by  DHS  for  the  fastest 
growing  mode  in  an  incompressible  flow.  It  turns  out  that  the  effects  of  compressibility 
may  be  scaled  out  with  the  result  that  the  eigenvalue  problem  determining  the  growth  rate 
of  the  fastest  growing  mode  is  identical  to  that  solved  by  DHS.  Thus  we  do  not  present 
the  system  of  equations  here  but  point  out  that  full  details  of  the  compressible  problem 
may  be  obtained  from  the  second  author.  Consequently,  the  results  of  DHS  also  describe 
the  solution  for  a  compressible  fluid.  We  have  that  each  unstable  mode  has  its  maximum 
growth  rate  occurring  at  a  finite  value  of  A*.  DHS  showed  that  the  most  unstable  mode 

“X  j2— — "^1 

corresponds  to  A  =  0.476,  ^  =  0.312  where  a  =  and  the  growth  rate  of  the 

fastest  growing  mode  is  given  by  G? The  results  described  above  show  that 
for  a  compressible  fluid,  as  well  as  an  incompressible  fluid,  the  most  unstable  linear  Gortler 
vortex  at  high  Gortler  numbers  is  viscous  with  wavenumber  of  O(G^),  rather  than  0(Gf), 
which  is  appropriate  to  the  unstable  modes  close  to  the  right-hand  branch  of  the  neutral 
curve.  An  important  result  is  that  the  most  unstable  modes  occur  close  to  the  wall.  This 
suggests  that  significant  coupling  coefficients  will  be  possible  in  the  receptivity  problem  for 
the  most  unstable  modes.  This  was  shown  to  be  the  case  for  the  incompressible  problem  by 
DHS  and  identical  rt  ults  for  a  compressible  fluid  may  be  inferred  simply  from  the  results  of 
DHS.  Of  interest  would  be  the  behavior  of  the  fastest  growing  mode  in  the  hypersonic  limit. 
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Figure  11.  The  position  rj  of  the  trapped  layer  modes  for  a  »  1  as  a  function  of  Moo  for 
adiabatic  fluid  and  also  for  an  isothermal  fluid  with  T^,  =  0.5, 1,2,4  8 
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Figure  12.  The  neutral  growth  rate  lio  from  (6.3)  for  the  trapped  layer  modes  as  a  function 
of  Moo  for  an  adiabatic  fluid  and  also  for  an  isothermal  fluid  with  =  0.5, 1,2, 4. 8. 
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Figure  13.  The  growth  rate  /?  as  a  function  of  a  for  =  5:  —  numerical  solution  of  (2.16); 
- asymptotic  solution  from  (6.6). 
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F'igure  14.  Sketch  of  the  near  crossing  of  the  wall  layer  and  adjustment  layer  modes.  The  con¬ 
tinuous  curves  represent  the  finite  Mach  number  situation  and  the  dashed  lines  the  situation 
in  the  asymptotic  limit  — +  oo. 
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